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Abstract 

A  simplified  open-loop  dynamic  model  of  the  T700  turboshaft  engine,  valid  within  the  normal  operating  range  of 
the  engine,  is  developed.  This  model  is  obtained  by  linking  linear  state  space  models  obtained  at  different  engine 
operating  points.  Each  linear  model  is  developed  from  a  detailed  nonlinear  engine  simulation  using  a  multivariable 
system  identification  and  realization  method.  The  simplified  model  may  be  used  with  a  model-based  real  time 
diagnostic  scheme  for  fault  detection  and  diagnostics,  as  well  as  for  open  loop  engine  dynamics  studies  and  closed 
loop  control  analysis  utilizing  a  user  generated  control  law. 

Nomenclature 

Variables 

A,  B,  C,  D,  E,  F,  G,  K,  Z 

N. 

Nmk 
Np 

Qr 
Qs 
T45 
u 

Wp 
x 

XCPC 

y 


state  space  matrices 
gas  generator  speed 
main  rotor  speed 
power  turbine  speed 
static  pressure  at  station  3 
requited  torque 
engine  shaft  torque 
interturbine  gas  temperature 
engine  model  input  vector 
fuel  flow 

engine  model  state  vector 

collective  pitch 

engine  model  output  vector 


Greek  symbols 
8[variable] 

It 


normalized  deviation  of  variable  from  nominal 
observability  index 


model 
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Introduction 


Full  nonlinear,  component-based  tuibine  engine  simulations  have,  generally,  been  large,  slow,  unwieldy  computer 
programs  which  requite  large  amounts  of  computing  power  and  time  to  produce  small  amounts  of  data.  Due  to  their 
size  and  complexity,  however,  these  nonlinear  simulations  are  difficult  to  use  for  various  control  design  and 
diagnostics  studies.  In  this  paper,  a  simplified  state  space  model  of  the  T700  engine  is  developed  which  can  be  used 
for  control  system  design  and  model-based  fault  detection  and  diagnostics  studies.  This  simplified  model  is 
developed  using  the  data  generated  from  a  component-based  real-time  simulation  of  a  T700  turboshaft  engine  which 
is  described  in  [1], 

Customarily,  the  development  of  simplified  real-time  turbine  engine  simulations  involves  creating  piecewise  linear 
state  space  perturbation  models  [2]  at  several  operating  points,  enough  to  cover  the  operating  range  of  interest.  The 
set  of  linear  systems  covering  the  operating  range  is  linked,  sometimes  using  a  curve-fitting  technique,  so  that  the 
resulting  equation  contains  matrices  whose  elements  are  described  by  polynomials  evaluated  at  the  current  operating 
point.  This  methodology  was  used  in  [3],  [4],  and  [5]. 

In  this  paper,  the  procedure  used  to  develop  a  simplified  dynamic  model  of  the  T700  turboshaft  engine  is  described. 
This  approach  was  successfully  applied  to  obtain  simplified  models  for  control  and  diagnostics  design  in  other 
propulsion  systems,  for  example  [6].  Linearized  dynamic  models  of  the  engine  at  five  different  operating  conditions 
were  obtained  by  using  a  multivariable  system  identification  and  realization  technique  [7]  with  the  input  and  output 
data  generated  from  the  nonlinear  simulation  described  in  [1].  The  simplified  model,  covering  the  whole  range  of 
operation  of  the  engine,  was  then  developed  by  combining  the  point  models.  The  accuracy  of  the  simplified  model 
was  checked  by  comparing  the  responses  of  certain  engine  variables  to  the  corresponding  responses  from  the 
nonlinear  simulation  for  an  input  signal  which  drove  the  simplified  model  through  a  range  of  operating  points. 

A  possible  application  of  this  model  is  a  model-based  diagnostics  scheme  for  fault  detection  purposes  which 
compares  the  parameters  of  the  simplified  model  with  the  identified  parameters  of  the  actual  system  on  a  continuous 
basis.  The  differences  between  the  modeled  and  identified  parameters  may  be  analyzed  to  detect  the  presence,  to 
isolate  the  source,  and  to  determine  the  extent  of  faults.  Other  applications  for  this  model  would  include  open  loop 
engine  dynamics  studies,  and  closed  loop  control  analysis  with  a  user-generated  control  law. 

The  next  section  of  this  paper  contains  a  description  of  the  T700  engine.  This  is  followed  by  sections  describing 
the  modeling  and  the  identification  of  the  open  loop  dynamics  of  the  engine.  Next,  results  obtained  from  the 
identified  linear  operating  point  models  are  compared  with  the  results  obtained  from  the  nonlinear  simulation  for  the 
same  input.  Finally,  simulated  output  from  the  simplified  model  is  compared  with  the  results  obtained  from  the 
nonlinear  simulation  for  a  given  test  input. 

Description  of  the  T700-GE-700  Engine  in  a  Helicopter  Application 

The  T700-GE-700  is  a  turboshaft  engine  used  in  pairs  in  the  Army’s  Apache  and  Blackhawk  helicopters.  This  1600 
horsepower-class,  modular,  two-spool  engine  (figure  1)  consists  of  a  gas  generator  section  and  a  free  power  tuibine 
[8].  The  gas  generator  section  is  made  up  of  a  five-stage  axial  and  a  single-stage  centrifugal  compressor,  a  low  fuel 
pressure  flowthrough  annular  combustion  chamber,  and  an  air-cooled,  two-stage,  axial-flow,  high-pressure  turbine. 
The  free  power  turbine  is  a  two-stage,  uncooled,  axial-flow  type.  There  exists  a  one-way  coupling  between  the 
power  turbine  and  the  gas  generator,  i.e.  the  power  turbine  extracts  work  from  the  gas  turbine  cycle  but  does  not 
otherwise  affect  it.  Through  mechanical  linkages  and  gears,  the  power  turbine  drives  the  rotor  system  of  the 
helicopter  such  that  the  rotor  angular  velocities  are  directly  proportional  to  power  turbine  speed.  Thus,  in  the 
helicopter  application,  it  seems  natural  to  consider  the  power  turbine  as  part  of  the  rotor  system,  to  which  it  is 
mechanically  linked  [9]. 

The  helicopter  pilot  commands  an  altitude  change  by  moving  the  collective  stick  which  alters  the  collective  pitch 
of  the  main  rotor  blades.  A  change-in  the  blade  pitch  angle  causes  a  change  in  load  as  lift  is  increased  or  decreased. 


2 


The  standard  7700  control  system  regulates  power  turbine  speed  with  fuel  flow  to  accommodate  variations  in  the 
load.  The  power  turbine  shaft  is  designed  to  operate  at  a  constant  20,900  rpm.  Changes  in  this  power  tuibine  speed 
are  induced  by  disturbances  such  as  collective  pitch  changes,  wind  gusts,  etc.  As  the  engine  reaches  a  steady  state 
condition,  the  control  system  causes  the  power  turbine  speed  to  return  to  its  design  point.  The  operating  point  of 
the  engine  is  defined  by  the  gas  generator  speed  as  a  percent  of  its  design  speed  as  in  [10]. 

As  seen  in  figure  2,  the  two  engines  work  in  tandem  employing  a  torque-sharing  arrangement  to  turn  the  rotor 
system.  In  each  engine,  the  power  turbine’s  shaft  is  connected  to  the  main  gearbox  through  a  nose  box  (a  small 
reducing  gearbox  which  redirects  the  torque  by  90  degrees)  and  a  freewheeling  clutch.  The  main  rotor  is  driven  by 
a  shaft  projecting  from  the  gearbox.  The  clutch  allows  the  turbine  to  drive  the  rotor  system  but  not  vice  versa.  TTie 
tail  rotor  is  driven  by  a  shaft  extending  from  the  transmission  which  is  turned  though  a  gear  mechanism  powered 
by  the  main  rotor.  Thus,  it  spins  as  long  as  the  main  rotor  is  spinning,  even  if  the  engine  is  disengaged.  When  the 
helicopter  tries  to  climb,  implying  that  torque  must  be  increased,  the  control  system  increases  the  fuel  flow  to  the 
engines.  The  rotor  applies  a  load  to  the  turbine  as  long  as  they  are  engaged.  A  command  to  descend  causes  the 
control  system  to  reduce  fuel  flow,  thereby  reducing  the  torque  applied  to  the  rotor  from  the  shaft.  The  rotor  will 
spin  freely  along  with  the  shaft,  but  will  only  be  driven  by  the  shaft  when  friction  slows  the  rotor  down.  This  leads 
to  a  situation  where  the  load  depends  on  the  direction  of  motion  of  the  helicopter.  A  command  to  climb  causes  a 
significant  load  and  the  associated  excursions  in  the  state  variables.  On  the  other  hand,  a  command  to  descend 
causes  a  sudden  drop  in  load  because  the  rotor’s  inertia  keeps  it  spinning  at  the  desired  velocity  even  though  fuel 
flow  is  reduced,  meaning  shaft  torque  is  decreased. 


Modeling  of  the  Open  Loop  Dynamics  of  the  Engine 

As  mentioned  earlier,  it  is  reasonable  to  consider  the  power  turbine  to  be  part  of  the  load  rather  than  the  engine 
because  of  the  one-way  coupling  and  the  direct  linkage  between  the  power  turbine  and  the  rotor.  Figure  3  shows 
a  block  diagram  of  the  open-loop  engine/rotor  system. 

The  basic  objective  of  the  control  system  of  the  engine  is  to  maintain  a  constant  power  turbine  speed  in  the  presence 
of  load  disturbances.  This  is  accomplished  by  modulating  the  fuel  flow  rate  and  the  variable  compressor  geometry. 
However,  for  the  purposes  of  this  study  and  consistent  with  the  nonlinear  simulation,  the  variable  geometry  is 
assumed  to  be  scheduled  open  loop  as  a  function  of  compressor  inlet  temperature  and  gas  generator  speed. 
Following  [9]  and  [11],  the  rotor  system  is  treated  as  a  separate  dynamical  system  with  engine  shaft  torque,  Qs, 
as  input  and  the  main  rotor  speed,  as  output  as  shown  in  figure  3.  The  main  rotor  speed  is  assumed  to  be 
proportional  to  the  turbine  speed  by  a  constant  gear  ratio  since  the  shaft  elements  between  the  hub  and  the  power 
turbine  are  nearly  rigid  [9], 


In  figure  3,  the  inputs  to  the  engine  are  the  power  turbine  speed,  NP,  and  the  fuel  flow  rate,  WF.  The  outputs  are 
the  gas  generator  speed,  Ng,  engine  torque  transmitted  by  the  power  turbine  shaft,  Qs,  compressor  static  discharge 
pressure,  PS3,  and  power  turbine  inlet  temperature,  T4J  (figure  1).  All  variables  are  normalized  by  dividing  them  by 
their  values  at  100%  design  gas  generator  speed.  The  pressure  and  the  temperature  are  included  in  the  output  since 
they  are  used  by  the  control  system.  Power  turbine  speed  is  considered  as  an  input  to  the  engine  in  order  to 
reference  the  engine  shaft  torque  with  respect  to  that  speed  and  to  take  into  account  the  frictional  losses. 


In  this  study,  the  identification  and  modeling  of  the  open  loop  dynamics  of  the  helicopter  engine  as  described  above 
is  considered  separately  from  the  rotor  system.  Once  the  model  of  the  engine  is  obtained  it  can  be  coupled  with  any 
desired  rotor  system  and  the  performance  of  the  combined  engine/rotor  system  can  be  analyzed. 


The  nonlinear  open  loop  dynamics  of  the  T700  helicopter  engine  can  be  described  by  the  nonlinear  differential 
equations 
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y(t)  =  g(x(t ))  (2) 

where  x,  u,  and  y  are  the  state,  the  control  and  the  output  vectors,  respectively.  Linearizing  these  equations  about 
a  nominal  operating  condition  and  discretizing  yields 


8x(*  +  1)  =  A  5x(k)  +  B  8u(k)  (3) 

by{k)  =  C  bx(k)  (4) 

where  8x,  8u,  and  8y  are  the  normalized  deviations  of  the  state,  the  input,  and  the  output  vectors  about  the  nominal 
operating  conditions  defined  as 
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where  the  subscripts  i,  n  and  0  refers  to  the  i*  component,  the  normalization  value  and  the  nominal  operating 
conditions.  The  normalization  values  and  the  nominal  operating  condition  values  are  given  in  Table  1. 


It  is  assumed  that  the  system  described  by  equations  (3)  and  (4)  is  stable  and  observable  and  the  C  matrix  has  full 
row  rank.  Furthermore,  this  system  is  realized  in  a-canonical  form  developed  in  [7]  and  [12],  i.e.,  the  following 
relations  hold: 


C  =  [0  :  H  '] 

A  -  Ao  +  KHC,  with  A*  =  0 

(6) 

woa:  =  o 

( HC)rA*Kc  =  0,  for  p.  >  p;  and  k  <  \i.  -  py 

Here  the  subscripts  r,  and  c,  denote  the  i*  row  and  j*  column  respectively.  Superscripts  indicate  exponentiation.  The 
structure  matrix  A,  is  lower  left  triangular  and  consists  of  zeros  and  ones  only  and  is  determined  by  the  observability 
indices,  p„  where  i  associates  p,  with  the  i*  output  and  p  =  max  { p, )  [13].  The  matrix  K  i.  a  deadbeat  observer  gain 
[14].  A  procedure  for  obtaining  a  state  transformation  matrix  to  represent  a  linear  system  in  a-canonical  form 
is  presented  in  the  appendix  together  with  an  example. 

System  Identification  Using  the  a-Canonical  Form 

System  identification  is  defined  [15]  as  the  determination,  on  the  basis  of  input  and  output,  of  a  system  within  a 
specified  class  of  systems,  to  which  the  system  under  test  is  equivalent.  In  this  study,  the  specified  class  of  systems 
considered  for  point  models  around  nominal  operating  conditions  is  the  class  of  linear  systems  described  by  equations 
(3)  and  (4). 

Selection  of  an  appropriate  input  signal  to  generate  the  output  data  is  an  important  step  in  identification  problems. 
A  basic  criterion  for  this  selection  is  that  the  input/output  data  should  be  informative  enough  to  discriminate  between 
different  models  among  the  class  of  models  being  considered  [16].  Without  this  discrimination  there  is  no  guarantee 
that  the  obtained  parameters  are  the  true  parameters  of  the  system.  This  criterion  can  be  expressed  mathematically 
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in  terms  of  the  covariance  matrices  of  the  input  signals  and  the  order  of  the  system  being  identified  [16].  Pseudo¬ 
random  multi-level  sequences,  such  as  those  used  in  this  study,  are  examples  of  input  signals  that  can  be  used  for 
designing  an  informative  experiment. 

Once  a  class  of  systems  and  the  input  signal  is  selected,  the  identification  problem  becomes  one  of  determining  the 
A,  B,  C  matrices  of  the  system  and  its  structure.  To  accomplish  this,  a  multivariable  system  identification  technique 
developed  in  [7]  based  on  the  a-canonical  form  is  used.  This  technique  is  briefly  explained  below  for  completeness. 

The  expression  for  the  state  §x(k)  at  time  k  of  the  system  described  by  equations  (3)  through  (6)  is  given  by 
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Using  the  nilpotency  of  A0  the  above  equation  yields 
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The  observability  indices,  pi;  are  evaluated  from  the  observability  matrix  constructed  from  the  input  and  output  data 
as  explained  in  [7],  If  the  observability  index,  p,  is  known,  then  equation  (9)  contains  the  unknown  A,  B  and  C 
matrices  for  given  inputs  and  outputs.  Using  a  data  length  of  N  and  a  least  squares  estimation  technique,  equation 
(9)  can  be  solved  for  the  unknown  parameters. 


Linearized  Point  Models  of  the  T700  Engine 


Using  the  above  modeling  and  identification  procedure,  the  parameters  of  the  linear  point  models  are  estimated  at 
five  different  operating  points  corresponding  to  different  percentages  of  gas  generator  design  speed  levels  as  shown 
in  table  1.  At  each  operating  point,  p  was  found  to  be  equal  to  2.  As  mentioned  earlier,  the  input  vector,  8u, 
represents  the  normalized  deviations  of  the  power  turbine  speed  and  the  fuel  flow  rate  from  the  nominal  values  (table 
1)  defined  as 


8u  =  [8NP  8W,]r  (10) 

The  output  vector,  8y,  represents  the  normalized  deviations  of  the  gas  generator  speed,  engine  shaft  torque, 
interturbine  gas  temperature,  and  the  compressor  exit  static  pressure  from  their  nominal  values.  It  is  defined  as 


Sy  =  [8JV  8 Qs  8T4,  8PJt  (11) 

The  validity  of  the  estimated  parameters  of  the  system  is  checked  by  comparing  the  responses  obtained  from  the 
identified  system  with  the  response  of  the  nonlinear  simulation.  The  model  of  the  identified  system  utilizes  the 
measurement  of  the  input  data,  8u,  to  predict  the  output,  8ym,  and  the  state,  8xm,  as 
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A  tri-level  pseudo-random  sequence  is  used  for  the  fuel  flow  input.  Since  the  power  turbine  speed  input  is  basically 
used  as  a  reference  speed  for  engine  shaft  torque,  a  signal  consisting  of  steps  and  ramps  which  cover  he  possible 
range  of  power  turbine  speed  is  used  for  that  input.  Both  of  the  input  signals  are  shown  in  figure  4. 

Most  of  the  interaction  between  the  rotor,  drive  train,  and  engine  takes  place  at  or  below  the  main  rotor  frequency 
of  about  300  rpm  (5  Hz).  Thus  it  is  important  that  the  simulation  be  accurate  in  the  range  below  5  Hz.  Based  on 
a  priori  knowledge  about  the  maximum  rise  time  of  the  engine,  a  clock  time  of  2  seconds  and  a  length  of  26  (the 
number  of  unique  outputs  of  the  pseudo-random  generator)  is  used  with  the  tri-level  sequence.  The  sampling  time 
is  selected  to  be  0.1  second.  This  corresponds  to  a  maximum  frequency  of  31.4  rad/sec  (5  Hz),  a  minimum 
frequency  of  0.12  rad/sec  (0.02  Hz)  and  a  signal  duration  of  52  sec. 

The  test  input  signal  consists  of  a  full  length  tri-level  pseudo-random  sequence  covering  the  same  frequency  range 
as  the  one  used  for  identification  as  shown  in  figure  5.  In  order  to  compare  the  responses  of  the  identified  system 
with  the  responses  obtained  from  the  nonlinear  simulation,  a  standard  error  of  estimate  (SEE)  is  defined  as 


E  [5y„,(*>  "  ty(*)f 

SEE  =  _  (14) 

£  byhk) 

\  bm\ 

Here  the  subscript  i  denotes  the  i81  element  of  the  output  vector  8y,  and  8ym  refers  to  the  model  output.  The  SEE’s 
obtained  for  the  test  signal,  given  in  table  2,  indicate  very  good  agreement.  Comparison  of  the  responses  of  the 
identified  model  to  the  responses  obtained  from  the  nonlinear  simulation  for  the  input  signals  shown  in  figure  5  at 
an  operating  condition  corresponding  to  96%  gas  generator  design  speed  are  shown  in  figures  6  through  9.  The  tests 
were  done  using  the  open  loop  nonlinear  and  linear  point  models  of  the  T700.  The  figures  show  how  closely  the 
linear  point  model  matches  the  nonlinear  simulation  for  perturbations  about  an  operating  point.  These  figures  are 
representative  of  the  responses  obtained  at  all  operating  points. 

A  Simplified  Model  of  the  T700  Helicopter  Engine 

A  simplified  model  of  the  T700  helicopter  engine  is  obtained  by  linking  the  linear  point  models  described  above 
through  look-up  tables  in  a  computer  program.  In  this  scheme,  the  gas  generator  speed  is  checked  after  each  time 
step  and  the  A,  B  and  C  matrices  corresponding  to  the  current  operating  conditions  are  used. 

The  computer  program  uses  the  actual  (not  normalized)  values  of  the  inputs  and  the  outputs  rather  than  the 
normalized  deviations.  Hence  the  five  point  models  used  are  defined  as 


x{k+ 1)  =  D  x{k)  +  E  u(k)  +  F  (15) 

y(k)  -  G  x(k)  +  Z  (16) 

The  values  of  these  matrices  for  each  of  the  operating  conditions  are  given  in  table  3.  The  order  of  the  system 
model  at  each  operating  point,  8,  equals  the  sum  of  the  state  variables  times  their  corresponding  observability  index 
at  that  operating  point. 
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The  linked  model  was  tested  for  accuracy  across  a  range  of  operating  points  by  comparing  its  response  with  that  of 
the  nonlinear  simulation  under  simulated  load  conditions.  A  simple  model  representing  the  dynamometer  described 
in  [11,  is  used  as  the  simulated  load.  The  load  is  variable,  based  on  a  simulated  collective  pitch  input,  XCPC.  The 
equations  defining  the  load  from  [1]  are 


np  -  ci  /«?,  -  <?*)  dt 


(17) 


Qr  =  C2  (C,  -  XCPC(C4  -  c5  XCPO)  N2p 


(18) 


where  Qs  and  QR  are  th-  engine  shaft  torque  and  the  required  torque.  The  values  of  the  constants  are:  C,=0.534753 
l/(lb-ft-sec2),  C2=1.86x10'7  lb-fit-sec2,  C,=75.0,  C4=0.1047,  and  C,=0.1085. 


To  test  the  validity  of  the  simplified  model,  the  profile  of  the  collective  pitch  input,  XCPC,  shown  in  figure  10  was 
used  as  input  to  the  dynamometer  model  (equations  (17)  and  (18)).  Beginning  with  the  initial  conditions  of  the 
engine  and  XCPC,  the  dynamometer  model  and  closed  loop  nonlinear  simulation  were  numerically  integrated 
simultaneously  to  produce  the  responses  represented  by  solid  lines  in  figures  12-16.  The  values  of  fuel  flow  rate, 
WF  (figure  11),  obtained  from  the  control  system  portion  of  the  nonlinear  simulation,  and  power  turbine  speed,  NP 
(figure  12),  obtained  from  the  nonlinear  simulation,  were  then  used  as  input  to  the  simplified  model.  Figure  12 
shows  the  comparison  of  the  NP  values  obtained  directly  from  the  nonlinear  model  and  from  the  simplified  model 
determined  by  using  the  torque  output,  Qs,  in  the  dynamometer  model  (equations  (17)  and  (18)).  As  shown  in 
figures  13-16,  the  simplified  model  passes  through  several  operating  regions  (as  defined  by  Ng  (figure  13))  and  the 
variables  match  those  of  the  nonlinear  simulation  closely. 

Conclusions 

A  multivariable  system  identification  technique  is  used  to  obtain  point  models  of  the  T700  helicopter  engine  from 
data  generated  by  a  nonlinear  simulation.  These  point  models  are  used  to  obtain  a  simplified  model  of  the  engine. 
The  validity  of  the  simplified  model  is  checked  by  comparing  the  response  of  the  nonlinear  simulation  with  the 
response  obtained  from  the  simplified  model  under  simulated  load  conditions.  The  comparison  indicates  good 
agreement  between  these  responses.  Therefore  it  has  potential  to  be  used  with  a  model-based  real  time  fault 
detection  and  diagnostics  design,  as  well  as  for  open  loop  engine  dynamics  studies.  Because  it  is  in  state  space  form, 
the  model  is  convenient  for  closed  loop  control  analysis  utilizing  a  user  generated  control  law. 

Appendix:  The  a-Canonical  Form 

Consider  a  linear,  time-invariant  system  with  an  irreducible  realization  given  as 


<7(*+l)  =  Dq(k)  +  Eu{k)  ri<y. 

y(k)  -  Fq(k) 

with  observability  indices  p,,  where  i  associates  p(  with  the  i*  output.  Here  the  state  vector,  q,  the  input  vector,  u, 
and  the  output  vector,  y,  are  nxl,  mxl,  and  pxl  respectively.  A  state  transformation  T  can  be  used  to  represent  this 
system  in  a-canonical  form  as 


A  -  IDT  1  ,  B  -  TE  ,  C  -  FT  1  ,  (20) 

where  the  A  and  C  matrices  satisfy  (6).  The  state  transformation  matrix  T  is  obtained  by  using  a  reduced 
observability  matrix,  V„  which  is  obtained  from  the  observability  matrix,  V.  after  deleting  all  rows  which  are  linearly 
dependent  on  the  previous  ones.  The  matrix  T  is  defined  [12]  as 


7 


( HF ) 


«.i 


Kt  A  £>//, 


(21) 


H{  A  [DMrlJ,  D"*'^  ....  D*>~ldJ 


H  A  p,)' ■'  J 

where  (HF), ,  is  the  matrix  obtained  from  HF  after  deleting  the  rows  with  associated  observability  indices  p,  <  t. 
The  column  vectors,  d,  denote  those  columns  of  V,'1  such  that 


FrDVr'dt  =  1  (22) 

rt  1 

where  the  subscript  r,  denotes  the  i*  row.  The  deadbeat  observer  gain,  K,  for  the  transformed  system  of  (6)  is 
obtained  from  that  of  the  untransformed  system,  K,,  by 

K  =  TKy  (23) 

Example:  Let  (D,EF)  be  an  irreducible  system  with  the  realization 


0 

1  0 

e, 

D  - 

*i 

*2  *3 

,  £  = 

«2 

* 4 

*5  *6 

*3 

1  0  0 
0  0  1 


(24) 


where,  from  V,  the  observability  indices  are  p,  =  2  and  p2  =  1.  For  this  system,  the  observability  matrix,  V,  the 
reduced  observability  matrix,  V„  and  the  vectors  d,  and  d;  are  obtained  as 


V  = 


1 

0 

0 


*1*5  +  *4*6 


0 

0 

1 

*j 
*2 
f  *2*5 


+  *5*6  *3*5  +*« 


K 


K.>  = 


1  0  0 
0  0  1 
0  1  0 


0 

0 

4 

1 

,  4  = 

0 

0 

1 

(25) 


(26) 


Hence  H,  K,  and  H,  can  be  evaluated  as 
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noting  that 


- 

1  0 

k2  0 

,K{  - 

k2  0 

A(>A2+A?A,  A3 

,  H  - 

1  0 

A,  1 

A,  1 

a4+a,a,+*,a«  a6 

I  1  0  0 

0 HF).2  -  II  0  01  ,  WF).X  -  \  .  n  . 

■  |  *,  0  1 

The  transformation  matrix  T  can  be  determined  as 


(27) 


(28) 


(29) 


(30) 
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Table  1.  Operating  Points  and  Normalization  Values  Used  for  Model  Development 


operating  point. 


%NP  (%)  100.0 

Wf  (lb/sec)  0.0884 
%N,  (%)  87.5 

Qs  (ft-lb)  0369. 

T45  (deg  R)  1472. 

PS3  (psi)  132.0 


100.0 

0.1170 

91.4 

15309. 

1578. 

161.7 


100.0 

0.1456 

94.5 

21443. 

1675. 

188.9 


100.0 

0.1742 

96.7 

27303. 

1778. 

212.3 


100.0 

0.2028 

98.5 

32619. 

1896. 

231.5 


normalization 


100.0 

0.217 

99.5 

35150. 

1954. 

240.5 


Gas  generator  design  speed  =  44,700  rpm 


Power  turbine  design  speed  =  20,900  rpm 


Table  2.  Standard  Error  of  Estimate  (SEE)  for  each  Point  Model 


operating  point:  1 


Table  3.  State  Space  Matrices  at  each  Operating  Point 


Operating 

Point 

2 


%N=91 


Operating 

Point 

3 


%Ns=94 


Operating 

Point 

4 


%Ng=96 


0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000 

1.0000  0.0000  0.0000 

0.0000  1.0000  0.0000 

0.0000  0.0000  1.0000 

0.0000  0.0000  0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 


1 


.0303 
.0263 
.0030 
0.0155 
0.0387  |  0.1179 
0.9875 
0.0224 
0.0027 


0.1121 

.3146 

.4131 

.1308 

.1333 

.3105 

.4496 

.1443 


G2/10,000 


Zj/1  0,000 


0.0000  0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000  0.0000  0.0000 

0.0000  0.0000  0.0000  0.0000  0.0000  I  0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 


G3/10,000 


1.2196 


ZV  10,000 


0.0095 

2.1443 

0.1675 

0.0189 


G^  10,000 


1.1848 


ZJ\ 0,000 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000  0.0000 
0.0000  0.0000 
0.0977  0.0000 


0.0212 


Operating 

D, 

I 

h 

F, 

Point 

5 

0.0000 

0.0000 

0.0000 

0.0000 

0.1907 

I 

-2.2727 

0.0000 

0.0000 

0.0000 

0.0000 

0.2123 

BBSS 

j 

mm 

-1.8996 

%Ng=99 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0201 

; 

D 

0.0000 

-2.7338 

> 

0.0000 

0.0000 

0.0000 

0.0000 

0.2010 

K 

-1.7053 

1.0000 

0.0000 

0.0000 

0.0000 

0.8479 

IW 

2.4826 

HE ! ' 

0.0000 

1.0000 

0.0000 

0.0000 

-0.1461 

2.1337 

0.0000 

0.0000 

1.0000 

0.0000 

-0.0120 

Hpl? 

2.8933 

|||||? 

0.0000 

0.0000 

0.0000 

1.0000 

-0.0855 

0.1515 

1.9276 

022 

G5/ 10,000 

ZJlOfiOO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0098 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.7575 

0.0000 

3.2619 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.1896 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

ED 

0.0241 

0.0232 

FIVE-STAGE  LOW  FUEL  PRESSURE 

AXIAL  FLOW  FUEL  FLOW  INPUT  FLOWTHROUGH  ANNULAR 

COMPRESSOR  \  COMBUSTION  CHAMBER 


SINGLE-STAGE  TWO-STAGE  AXIAL  TWO-STAGE 


OUTPUTSHAFT  CENTRIFUGAL  FLOW  HIGH-PRESSURE  AXIAL  FLOW 

COMPRESSOR  TURBINE  POWER  TURBINE 

Figure  1 . — Cross  section  of  a  T700  turboshaft  engine. 
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Figure  2. — Twin  engine/rotor  system  in  a  helicopter. 


XCPC= collective  pitch 


Figure  3. — Open-loop  engine/rotor  system. 


-  oi  power  turbine  speed,  Np 

.  u 2  fuel  ^ow  rote,  Wp 


Time  (seconds) 


Figure  4. — Input  signals  u-j  (Np)  and  u2  (WF). 


'  -  N<j  from  nonlinear  model 

.  Ng  from  linear  point  mode 


Qs  from  nonlinear  model 
Qs  from  linear  point  model 


Figure  7. — Output  torque  at  96  percent  gas  generator  design  speed. 
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Power  turbine  speed  (7.  design  speed)  W r  (Ib/sec)  XCPC  (Percent  of  mcximun  volue) 


50 


Figure  1 0. — Collective  pitch  input. 


-  7.  Np  from  nonlinear  model 
7,  Np  from  simplified  model 


Figure  12. — Comparison  of  nonlinear  model  with  simplified  model. 
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Figure  15. — Comparison  of  nonlinear  model  with  simplified  model. 


Ps3  from  nonlinear  model 
Ps3  from  simplified  model 


Figure  16. — Comparison  of  nonlinear  model  with  simplified  model. 
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